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1 Introduction 



Interplay between statistical mechanics and combinatorics is one of the most interesting 
features in mathematical physics pp. Recently, a connection between the integrable phase 
model and enumeration of plane partition was shown in [Oj. The integrable phase model 
is a special limit of the g-boson model which is also related to the XXZ model fSJ ^1 El , 
and it can be solved applying the Quantum Inverse Scattering Method (QISM) [Tn| l9]. 

Study of generating function of plane partition |S1 IH E] has attracted many researchers 
going back to P. A. MacMahon. For an illuminating story of plane partition involving 
Alternating Sign Matrix theorems and the six-vertex model with domain wall boundary 
condition ^2 E] ? we refer to j2] . Recently, correlation functions for random (skew) plane 
partition were studied in a more general setting in terms of the free-fermion fields [3 
IH]. There the generating function of unrestricted plane partition is presented as just a 
normalization constant, or a zero-point function. 

In this article, we generalize [Oj and argue the connection between the integrable phase 
model and the boxed skew plane partition. The generating functions of the boxed skew 
plane partition are given by a special parameter choice of the generalized scalar products 
of the phase model. We calculate the generalized scalar products with generic parameters 
in two ways. Namely, the one is expressed in terms of the Schur function; the eigenfunction 
for the phase model, and another is, through the relations from QISM, expressed in terms 
of a determinant of matrix with elements of two-point functions and one-point functions. 

The paper is organized as follows. In Section |2l we introduce the phase model as a 
quantum integrable system. The integrability of this model is understood in the context 
of the QISM. Through graphical interpretations, the action of the monodromy operators 
on a Fock state turns out to be written in terms of the skew Schur functions. A scalar 
product, which is defined as the inner product of the A^-particle states of the phase model, 
is generalized. The inhomogeneous phase model is also considered. In Section IHl we briefly 
introduce the skew plane partition in a box and its generating function. In Section HI we 
see the connection between a generahzed scalar product of the phase model and the skew 
plane partition. A sequence of operators in a generalized scalar product corresponds to the 
shape of the skew plane partition, and the number of lattice sites is equal to the height. The 
special choice of the spectral parameters naturally gives the generating function of the skew 
plane partition through graphical observations. A generalized scalar product is calculated 
in two ways in Section El and Sectional The first one is by the help of the expression in the 
skew Schur functions. Another one is an inductive application of the commutation relations 
of the monodromy operators. In Section [3 we show that the generating function of the 
boxed plane partition has a determinant expression, which is a by-product of the explicit 
expression of a generalized scalar product obtained in Section (HI Section |H1 is devoted to 
the conclusion. 
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2 The Integrable Phase Model 



In this section, we introduce the integrable phase model and review the results on this 
model concerning its relation to the skew Schur functions based on [H]. We also introduce 
and generalize the scalar product of the wave functions of the phase model. This is the 
main quantity which is viewed as the generating function of the plane partition in the 
subsequent sections. 

2.1 Integrability of phase model 

We introduce the phase model as an integrable system. Consider a one-dimensional lattice 
of length M + 1 and label each site hj i = 0, ■ ■ ■ ,M. Let (pi, (pi denote the operators 
satisfying the following commutation relations 

[Ni,(pj] = -(piSij, [Ni, (p]] = (pl6ij, [(pi,(p]] = iTi6ij (1) 

for i,j = 0, ■ ■ ■ , M where Ni is the number operator and vTj is the vacuum projector. The 
L-operator for the phase model is defined as 

^•(") - ( i ) p) 

with the spectral parameter n G M. The local Fock states at site j are created from the 
local vacuum state |0)j such as \nj) = {(p])"'^\0)j,Nj\nj)j = nj\nj)j. The states satisfy 

0i|O)i = O, 

(l)j\nj)j = \nj-l)j, (l)]\nj)j = \nj + l)j, (3) 

The inner product of two Fock states is normalized, i{m\n)j = SijSnm- 
Let us introduce the i?-matrix (a 4 x 4 matrix) of the form 

Ru{u,v) = Ru{u,v) = f{v,u), 

i?22(M, v) = i?33(M, v) = g{v, m), (4) 

R23{U,V) = 1, 

where 

f{v,u) = —^ — 2^ ai^^u) = ^ , (5) 

U — V U — 

and the other elements are zero. Then we can find that the L-operator satisfies the inter- 
twining relation 

R{u,v){Ln{u) ® Ln{v)) = {Ln{v) ® Ln{u))R{u,v) . (6) 
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The monodromy matrix is defined as 

T(n) = LmLm-1 ■■■Lo=(^ ^ J^j ^ . (7) 

Hereafter, we refer to the elements of the monodromy matrix as the monodromy operators. 
The commutation relations of the monodromy operators can be obtained from 

R{u, v){T{u) ® T{v)) = {T{v) ® T{u))R{u, v). (8) 

Several explicit relations are listed in the following. 

[Aiu), Aiv)] = [Biu), Biv)] = [C{u), C{v)] = [D{u), D{v)] = 0, (9) 
giv, u)Ciu)Aiv) = fiv, u)Civ)Aiu), (10) 
g{v, u)D{u)B{v) = f{v, u)D{v)B{u), (11) 
C{u)B{y) = g{u,v){A{u)D{v) - A{v)D{u)} = g{u,v){D{u)A{v) - D{v)A{u)}ll2) 
C{u)A{v) = f{v, u)A{v)C{u) + g{u, v)A{u)C{v), (13) 
D{u)B{v) = f{v, u)B{y)D{u) + g{u, v)B{u)D{v), (14) 

Note that these commutation relations can be reduced from those of the g-boson model. 
2.2 Graphical representation of operators 

The monodromy operators A, B, C and D are expressed as the linear combination of prod- 
ucts of elements of L-operators. For later convenience, we introduce a graphical represen- 
tation of the L-operator and the monodromy matrix in this subsection. 

We assign a vertex with two vertical arrowed edges to each element of L-operators. 
Each arrow is directed upward or downward. The vertex is indexed by i for the i-th site. 
Pairs of arrows stand for (pi, u and respectively (see FigureH]). 

This vertex representation is useful for a calculation of products of L-operators. We 
attach vertices vertically by the edge when the direction of the arrow between two neighbor- 
ing vertices are the same. This gives a one-dimensional vertical graph. Then, each element 
of a product of L-operators can be represented as summation over graphs of all arrow con- 
figurations with certain fixed boundaries. For instance, the (1, l)-element of Li{u)Lj{u) is 
-|- and the corresponding graphs have both boundary arrows upward. 

The monodromy operators A, B, C and D have also one- dimensional graph representa- 
tions. From their definitions ((Tj), each operator can be represented as a sum over graphs 
of all possible arrow configurations on a vertical lattice with M + 1 vertices and fixed 
boundary arrows. Here the lattice sites are numbered by z = 0, ■ ■ ■ , M from the bottom. 
Operator B{u) (resp. C{u)) has the top and bottom arrows pointing outward (resp. in- 
ward). Operator A{u) (resp. D{u)) has the top and bottom arrows pointing upward (resp. 
downward) . 



4 



4 



(f)j^ u u 



Figure 1: Vertex representation for the elements of L-operator. 

The number of 0^ is bigger (resp. less) by one than that of in every possible arrow 
configuration of B{u) (resp. C{u)). In this sense, the operators B{u) and C{u) have the 
property of the creation and annihilation operator, respectively. On the other hand, the 
number of 0^ is exactly equal to that of in every possible arrow configurations of A{u) 
and D{u). 

In later discussion, we call those shown in Figure El the basic configurations. The other 
arrow configurations can be obtained by flipping some arrows upside-down. 

M 

M - 1 
1 


A B C D 

Figure 2: Basic arrow configurations of the operators A, B, C and D. 

2.3 Fock states and partitions 

In this subsection, we see the actions of the monodromy operators A, 5, C and D on 
a Fock state. We introduce a bijection from an arrow configuration considered in the 
previous subsection to a lattice path configuration. By using this bijection, we express 




each matrix element of the operators with respect to the Fock states in terms of the skew 
Schur functions. 

Any state of the phase model can be expanded by a linear combination of states of the 
form 

M M 

= JJ |n,),-, N = Uj > 0, (15) 

i=o j=o 

where {rij} stands for a configuration {rij} = {uq, ■ ■ ■ , um) with rij > 0. We call this kind 
of state an iV-particle basis. In particular, we denote the vacuum state by |0) = YljLo 
As the action of the operator B{u) (resp. C{u)) on a state increases (resp. decreases) 
the total number of particles by one, we call the state H^i -^("^jOlO) iV-particle state 
and {0\Ylj=iC{uj) its conjugate. Obviously the A^-particle state is a linear combination 
of A^-particle bases. A matrix element {{n'j}\B{u)\{nj}) vanishes unless \n'j — nj\ < 1 and 
A^' = A^ + 1. Similarly, a matrix element {{n'j}\C {u)\{nj}) vanishes unless \n'j — nj\ < 1 and 
A^ =A^-1. Since the actions of A{w) and D{w) preserve the number of particles, matrix 
elements {{n'j}\A{w)\{nj}) , {{n'j}\D{w)\{nj}) vanish unless \n'j — Uj I < 1 and A^ = A^. 

Now, we rephrase the arrow configuration of vertical graphs in the previous section 
as a non-intersecting up-right path configuration over the lattices. The action of the 
monodromy operator on a Fock state can be represented as a lattice path configuration on 
a vertical lattice. 

Let us fix an A^-particle basis {n'j} and an arrow configuration of a vertical graph. Put 
n'j horizontal paths entering the j-th site from the left of the lattice. Put also another 
path entering the 0-th site from the bottom if the bottom arrow is inward. The rule is 
that only one path can go along an upward arrow and no path along a downward arrow. 
If the j-th site and the (j + l)-th site are connected by an upward arrow, one of the paths 
goes upward from the j-th site along the vertical lattice until it faces a downward arrow 
or another horizontal path, then it goes rightward. Other paths getting into the j-th site 
just pass rightward horizontally. A path entering the M-th site and going up along the top 
arrow, when the top arrow is upward, just goes vertically. Note that if there is no vertical 
path on a certain upward arrow, every is annihilated by one of 0j's. As a result, we 

have a configuration {nj} by counting nj paths coming right out of the j-th site. When the 
monodromy operators act successively on a Fock state, we can construct up-right paths 
by repeating the above procedure. In this way, a lattice path configuration corresponds 
one-to-one to an arrow configuration of vertical graphs. Figure El is an example of this 
correspondence. 

We rewrite a Fock state in terms of a partition in order to obtain an explicit expression 
of the matrix elements of the monodromy operators. For a given A^-particle basis we 
assume that rij^ > {r = 1, ■ ■ ■ , R) and M > ji > j2 > ■ ■ ■ > Jr > 0. The identification is 
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Figure 3: An example of a lattice path configuration corresponding to an arrow configura- 
tion. 



that A = (Ai, ■ ■ ■ , Atv), M > Ai > A2 ■ ■ ■ > A^v > where 

Afc = jl, 1 < A; < Uj^ 

Afc = j2, rij^ + 1 < k < rij^ + rij^ 



R-l R 



(16) 

Afc = JR, + 1 < < ^ rij^ = N. 

r=l r=l 

We write the A^-particle basis as |A). The Fock space is spanned by the orthonormal 

basis I A) of partitions. In a lattice path configuration of the monodromy operator, the site 
Aj is exactly the passing site of the z-th path from above. 

Then, the matrix elements of monodromy operators with respect to |A) bases have the 
following expression in terms of the skew Schur functions Sa/^^: 

Proposition 2.1. Matrix elements of the monodromy operators are given as follows (A C 

{fi\A{u)\X) = u~^''+'^s,/,{u'), for fi C M^, 

(/i|5H|A) = u-''s,^,iu'), for /i C M^+\ 

{fi\Ciu)\X) = u^'s,/^,{u~'), for/iCM^-\ 

(^|D(n)|A) = u''+h,/,{u-'), for /i C M^. 

Proof. We consider a non- vanishing {fi\C{u)\X). In the view of a lattice path, the partitions 
/i and A are identified with the initial and final positions of the up-right paths, respectively. 
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Fix fi C M^~^ and consider the corresponding path configuration with respect to A C 

. We treat A C because another path comes in along the bottom arrow for the 
configuration of C{u). 

When, graphically, the arrow configuration for C{u) is the basic one, the path config- 
uration is such that all — 1 paths from the left pass through without going vertically, 
and another path comes in from the bottom and goes out at site 0. Then, A is equal to /i. 
Note that Aat = is irrelevant to the shape A. We find in this case {fi\C{u)\X) = u*^. 

For the case where X ^ fi, the arrow configuration is the one where some downarrows 
except the boundary ones are fiipped upside-down from the basic configuration. If one dow- 
narrow is fiipped from the basic configuration, we have another factor Let us consider 
a lattice path configuration where the i-th path (from above) goes along /j > uparrows. 
There are I = Yl^i uparrows fiipped upside-down from the basic configuration, and we find 
{fi\C{u)\X) = Now, the partitions /x and A satisfy the following properties. First, 

A 3 /i and Aj = /ij + U. Second, the skew partition A//i is a one-vertical strip since the 
condition that only one path should pass through every uparrow tells Aj+i < /ij < A,, Wi. 
From the definition of the skew Schur function, sa/^(m~^) = Together with these 

observations, we finally obtain (/i|C(n)|A) = u'^^ sx/fi{u~'^). 

Just similar considerations give matrix elements for the other operators. □ 

By inserting the complete set 1 = between every successive operators and 

using J2u ^^^/'^(^)^v/x{y) = s^/x{x,y), we furthermore have the following relations. 

Proposition 2.2. The matrix elements of products of the monodromy operators are as 
follows (X,X' CM^^fiC M^+-^). 



{fi\Y[B{uj)\X) = (mi---mj) ^^s^/a(mi,--- ,Mj), 

i=i 
J 

(A'l n ^(^.)|A) = (ni ■ ■ ■ u,)-(^^+i).y/,(n?, ■■■,u% 

(A| n ^K)|A') = («i ■ ■■ujr+h,,/,{u^\ ■ ■ ■ , uf). 



Other elements are zero. 

These expressions allow us to write a generalized scalar product as a sum of products 
of the skew Schur functions (see Section E)). 

The hermitian conjugates of the monodromy operators satisfy the following. 
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Proposition 2.3. The hermitian conjugates of the operators are related as 

A\u) = Din-'), D\u) = A{u~'), 

B\u) = C{u~^), C\u) = B{u~^), ^ ' 

i.e. for the monodromy matrix it holds that 

T\u) = T{u-^). (20) 



Proof From jllj), we have (/i|A(n)|A) = n-(^+^)s^/A(n2) = {\\D{u-^)\^i) and (/i|5(n)|A) 



u 



-^VA(^') = (A|C(n-i)|^). □ 



2.4 Generalized scalar product 



We briefly review tlie scalar product of tlie pliase model. The state IljLi is the 
A^-particle state. The inner product of A^-particle states is called a scalar product of the 
phase model. A scalar product of the phase model was first calculated in through the 
partition function of the six- vertex model with domain wall boundary condition |T!H IT^ . 
Another way to calculate it is to use the Schur function representation of the A^-particle 
state. The explicit expression of the scalar product of the iV-particle state is jS] 

M 



N N / N 

Ui 



{^\\{c{v,)\{b{um = 5^ .A(«?,---,OsA(^r 



J=l 3=1 \i=l / XCMN 



2 -2 

' ■ ■ ■ y'^N 



n Uj) Yl gi^i, Vj) det [H{vi, Uj)] i<^j<^, (21) 

i<j i>j 



where 



M+N / X M+N ' 

E[v,,u,) = {\^\ -(^1 (22) 




To prove (j2T|) . we just use the definition of the Schur function s a (x) = nj>j(^j~^«)~^ ^et x 
when the length of A is equal to or less than and the Cauchy-Binet formula. 

We define a generalized scalar product by the vacuum expectation value (VEV) of a 
sequence of the monodromy operators. It has a non-zero value only when the numbers of 
B and C are equal. In this article, we only consider the generalized scalar product of a 
Ci?-type word (see Section UTTjl : 

(oi n c(«f") n ^("^) ■ ■ ■ n n ^(^^) 

i=l i=l i=l i=l 

n ^(«?^) n ^("f ■ ■ ■ n ^("^) n ^iuf-m (23) 

i=l i=l i=l i=l 



N~j 
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A generalized scalar product of this type can be calculated in two ways: by use of the skew 
Schur function representation of the phase model, Proposition^! (see Section E}, and by 
use of the commutation relations of the monodromy operators (see Sectional). In both 
cases, a generalized scalar product can be expressed in a determinant form. 

A scalar product of the phase model is related with the boxed plane partition by 
adjusting the parameters as shown in jB]. Below, we show that a generalized scalar product 
is related with the boxed skew plane partition. A sequence of operators appeared in 
determines the base shape of the boxed skew plane partition. These subjects are considered 
in the successive sections. 

2.5 Inhomogeneous phase model 

We can generalize the phase model to that with inhomogeneous parameters; the spectral 
parameter u of Lj is replaced by uzj in each lattice site j. Since f{vzj, uzj) = f{v, u) and 
g{vzj,uzj) = g{v,u), the i?-matrix is unchanged so that the relations for the monodromy 
matrix obtained from the QISM still remain the same. The only change that comes into 
our game is the functions, 

M M 

a{u- {.,}) = n ^7'' d{u- {z,}) = u''^' n (24) 

The calculation procedure for the generalized scalar product in this inhomogeneous case 
is much easier through the QISM than the explicit eigenf unctions in terms of the Schur 
function. Only we have to do is a replacement ()24j) in the final expressions obtained for 
the homogeneous case. 

3 Plane Partition 

3.1 Boxed plane partition 

A plane partition is an array of non-negative numbers 

TT={7rij), 2,j = l,2,---, 

that is non-increasing in both i and j. A plane partition is sometimes called a 3D partition. 
Piling TTij cubes over the {i,j) square in the two-dimensional square lattice gives a three- 
dimensional object. The volume of vr is denoted by |7r| = '^^^^ ^ ^ boxed plane 
partition (or a plane partition in a box) iil<i<N,l<j<L and |7rjj| < M for some 
A^, L and M. A boxed plane partition can be considered as stacks of cubes in a box with 
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side lengths A^, L, M. The generating function of the plane partition in an x L x M-box 
is known to be 

N L ^ i+j+M— 1 

z,(iv,L,M)= g'^'=nn r' > (25) 

TveB{N,L,M) i=l 3=1 ^ 

where B{N, L, M) stands for the N x L x M box. If we take the limit N, L, M oo with 
the assumption |g| < 1, the partition function Zg converges to the celebrated MacMahon 
function, nr=i(l 

As explained in Section \2M\ a scalar product of the phase model can be graphically 
interpreted in terms of lattice path configurations. A lattice path configuration is also 
viewed as a 2D-projection of a boxed (skew) plane partition through the lozenge tiling (see 
Section 1^^ : a plane partition in an A^ x L x M-box is also equivalent to a lozenge tiling 
of an (A^, L, M)-semiregular hexagon. Hence there is a one-to-one correspondence between 
a boxed plane partition and a lattice path configuration. 

The generating function of the boxed skew plane partitions can be naturally obtained 
by a special choice of parameters of the generalized scalar product. Up to some factor, 
(0| IljLi ^("^j) rifcLi B{uk)\0) with Uj = , Vj = q~i, 1 < j < N becomes the generating 
function of the plane partition in an A^ x A^ x M-box. Taking the corresponding generalized 
scalar product, an explicit determinant expression of is given in Section [7| 

3.2 Boxed skew plane partition 

We introduce a skew shape X/fi where /i C A. A skew plane partition vr whose shape is A//i 
is an array of non-negative integers (vrjj) defined on the (2,j') -squares contained in a skew 
shape X/fi that is non- increasing in both i and j. The volume of a skew plane partition is 
|7r| = ^TTjj. We call n a boxed skew plane partition when A,/i C A^^ and \7Tij\ < M for 
some A^, L and M. 

4 Skew Plane Partition and Generalized Scalar Prod- 
uct 

4.1 Skew shape from a sequence of operators 

The scalar product of the A^-particle states is relevant to the plane partition in an A^ x 
A^ X M-box. The base of the box is an A^ x A^ square, corresponding to a sequence of 
the operators Y[f=i ^h^j) Y\^=i B{uj). More generally, a sequence of the operators for the 
generalized scalar product reads the shape of the boxed skew plane partition, as we will 
explain below. 
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Let us introduce a word consisting of four letters A, 5, C and D. This is what ignores 
parameters of A{u), ■ ■ ■ ,D{u) and vectors (0| and |0) in the generahzed scalar product. 
Let us define a C-B-type word by a word in which the first letters are A or C, and the 
next N' letters are B or D {N and A^' are non-negative integers). CACADDBB is an 
example of a CS-type word. 

A generalized scalar product of a C-B-type word corresponds to the plane partition 
with a skew shape A//i where A is a rectangle. We consider a CS-type word in which the 
numbers of A, B, C and D are A^^^, A^, A^ and N^, respectively. Note that the numbers of B 
and C are equal to have a non-zero value for the generalized scalar product. The rectangle 
A turns out to be (A^ + A'^i) x (A^ + A^^) (see Figure IH). We put this rectangle on the two- 
dimensional lattice in a way that two diagonal corners are (0, 0) and (A^ -|- A^^, A^ -|- Nd). 
Then, we draw a path as follows: 1) a path starts from the origin, 2) we read the word 
from left to right, and move one unit upward (resp. rightward) if the letter is either C or 
D (resp. either A or B). The shape which is surrounded by two edges of the rectangle 
and the zig-zag path and contains the point (0, A^ -|- Njj) determines fi. This is the one- 
to-one correspondence between a C-B-type word and a boxed skew shape. Explicitly, for a 
generalized scalar product of a C-B-type word the partition is given by 



--N + Na-'^ Km+i-p, 1 + XI -^"i+i-p - ^ - Lm+l-p, l< s <m, 



p=l 



p=l 



p=l 



(26) 



=A^ 



Na- K -J^Mg, L + 1 + X A^, < J < L + X A^g, 1 < r < n. 



9=1 



q=l 



5=1 



where we set N 



(0, N + A^^) 



YZi Li and K = YZi K^ = N. 

^ {N + Na,N + Nd) 







-^m—l 










Lm.—\ 




Mn-l 


Nn-2 






Mn 


Nn-l 






Nn 











(0,0) 



{N + A^^, 0) 



Figure 4: The skew shape determined by a sequence of the monodromy operators. 



12 



4.2 Lozenge tiling and skew plane partition 

In this subsection, we consider a map from an arrow configuration on vertical graphs in 
Section 12.31 to a lozenge tiling of a hexagon. A generalized scalar product with a special 
choice of the parameters can be naturally identified as the generating function of the skew 
plane partition. 

Let us consider a lattice path configuration of a CB-tjpe word. We arrange a lozenge 
tiling from this lattice path configuration. The domain of tiles is determined by the shape 
of the corresponding skew plane partition. Let one unit of a path going upward correspond 
to a tile of type a, and let one unit of a path going rightward correspond to a tile of type 
b in Figure The remained spaces are filled with tiles of type c. Figure IHl is an example 




a b c 



Figure 5: Three types of tiles. 

of the correspondence among a lattice path configuration, an arrow configuration and a 
lozenge tiling. We also draw paths on the lozenge tiling to help understanding. These 
paths on a lozenge tiling are indeed non-intersecting by construction. We can see that only 
non-intersecting lattice path configurations are mapped to lozenge tilings. 

A lozenge tiling gives a 2i5-projection of a 3D piling of cubes in a box, therefore, a 
boxed skew plane partition. 

In the arrow graph representation, we can show that the generalized scalar product of 
a C-B-type word is written as 

Ma Lb Ka 

E n iK-n^)"""" n n i^^r-''-^'^ n w^'f^-" n n(-"^)"^"^""^ (27) 

c n=l a m=l b 1=1 d k=l 

where the summation is taken over all arrow configurations, and ^ (resp. counts the 
number of up (resp. down) arrows except the boundary arrows in the i-th vertical graph. 
The powers of parameters are obtained by flipping some arrows upside-down from the basic 
configurations. Since the arrows at the boundaries are different, the powers of u'" and 
(or and u^) are different by one. 

Suppose that A is an (A^ + N^) x {N + Na) rectangle, is a partition (PBj) and vr 
is a skew plane partition whose shape is X/fi. Using the correspondence of a skew plane 
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Figure 6: An example of identification of a skew plane partition with a lattice path con- 
figuration and a lozenge tiling. 
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partition and an arrow configuration, we obtain 

N+Na N+Nd 



TT 



j=i j=i 



For 1 < ? < n, let s{Ai) be the set of Mi successive integers from X]j=i ^i + Sj=i + 
= U^^is{Ai), and let s(Cj) be the set of Ni successive integers from Mj- + 



E5=1 + 1' ^(C) = ur=is(Ci). Note that U s{C) = {1, 2, ■ ■ ■ , iV + A^^}. This de- 



composition of integers corresponds to the product chain Y[ " CY[ "^■■'11 
Similarly, for 1 < i < m, let s{Di) be the set of Li successive integers from X]j=i -^j + 
Y^^~}iKj + 1, s{D) = U^is{Di), and let s{Bi) be the set of Ki successive integers from 
Ej=i Lj + E5=1 ^^j + 1, = UZis{Bi). Note the relations. 



E E (j-i) = iA//ii, 

E ^' + E - 1) = l/^l + ^(^-4 - ^^)(^^ -Nd + 1] 
Then from ()28|) . we can show the following relation. 

5^ (-|)(M - 2ti) + 5^ (-|)(M - 1 - 2t|) 



(29) 



- E C-^){n-M)+ E C-^){2ti-M + r 



= kl + ^^^l/^l - yl^l + ^^(^^ - ^dKNa -Nd + 1). (30) 
We parameterize u^,u^, u'-^ and as 

{m^'} = {g"2|j e s{Ai)}, ioT 1 <i <n 



{u^'} = G s(-Bj)}, for 1 < z < m 

{w*^'} = {<l~^j e ■s(Ci)}, for 1 < i < n 
{n^^*} = {q^\j e s{Di)}, for 1 < z < m. 



(31) 



We denote by Zcspio*) the generalized scalar product with the parameterization (jHH). 
Combining (jTTj) and (jHUj) , we find the generating function of a boxed skew plane partition 
is equal to 

g-ZG5p(g-), 
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where 



OL = 



2M + 1 
2 




M + 1 
4 



(iV^-iV,5)(iVA-iVD + l). 



We give a remark about the upside-down property. The hermitian conjugate of a 
generahzed scalar product of a Ci?-type word yields the generating function of another 
kind of plane partition whose base shape is a hook Young diagram. From Proposition 12.31 
the hermitian conjugate of the monodromy operators make the arrow graphs upside-down 
and make a lattice path configuration and a lozenge tiling rotate 180 degrees. See Figure 
El rotated 180 degrees. 

5 Genelarized Scalar Product via Skew Schur func- 
tions 

5.1 Cauchy-Binet formula 

In the subsequent, we present a formula for the generalized scalar product that is a de- 
terminant of the skew Schur functions. For this purpose, the basic is the Cauchy-Binet 
formula for calculating determinant. 

Let two matrices A and B be of size m x n and n x m, respectively (without loss 
of generality, n > m). Let Smn denote the set of strictly increasing sequences of length 
m that can be chosen from {1,2,-- - ,n}. For any matrix X of size n x m and any S = 
{si, ■ ■ ■ , Sm} G Smn, we deuotc by X{S\m) the mxm matrix obtained from X by choosing 
all m columns of X and the m rows numbered by S, i.e. {X{S\m))ij = Xs^j. Similarly, if 
X is m X n, denote the mxm matrix Xi^.^l <i,j < m) hj X{m\S). The Cauchy-Binet 
formula shows that the determinant of the product of two rectangular matrices A and B 
is 



5.2 Skew Schur functions 

For any partition A = (Ai, ■ ■ ■ ,AAr),Ai > and some non-negative integer L, we define 




(32) 



A = (Ai, ■ ■ ■ , Xn+l) as 
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For any partitions A = (Ai, ■ ■ ■ , Atv) and /i C A, the skew Schur function has a determinant 
expression of the form [Jacobi-Trudi identity) 



det (/^A,^^,-^+i)l<i,■<^ 



det h 



:i,j<N+L 



(33) 



where h^, n G Z is the complete symmetric function. Applying the Cauchy-Binet formula 
!2|l. we can show the following expressions. 



Proposition 5.1. For non-negative integers N, L and M, we have 

'N+M+L 



^ sx/f,{x)sx/,y{y) = det 
^ s^,/xix)s^/x{y) = det 



hi_/iXx)hi_i^^{y) 



. 1=1 
'N+M+L 



l<i,j<N+L 



1=1 



(34) 
(35) 



l<i,i<N+L 



More generally, let A* = (A^^, ■ ■ ■ , AJy.), < ? < 2n and let = {x{, ■ ■ ■ , a^l^.), 1 < j < 2n 
denote 2n sets of variables, then we have a summation of products of the skew Schur 
functions over a = {X^\X\ < M, l<i<2n— 1} as a determinant of complete symmetric 
functions. 

n-l 

n ■5A2i+i/A2i(a;^*"'"^)sA2i+i/A2»+2(a;^*"^^) = det A-^ox-2^i{x^ , ■ 
where N = max{A^j|l <i< 2n},L = 'Y^Li, and 

N+L+M n-l 



,X 



2n\ 



l<i,i<N+L 



(36) 



2k+2\ 



(37) 



'i,--- ,hn-i=l k=0 

Let us take the limit M — 00. The generating function of Aij{x) becomes 



Y z'w^Mx\ 



2n\ 



i,j=l 



2n 



XjZ 



i=l j 



' n nti 

i=n+l j 



(3J 



where we take only positive powers and a{i) = (—1)*"^^. Here, suppose that two functions 
(1 — ^C^"*^)"^! C = z,w represent the formal power series of (: (1 — ^C^^)^^ = Xl^o '^"'C^"'- 



17 



If we also take partitions A , A = 0, we have 

n-l 



JJ^ SA2i+i/A2i(x^*'''^)sA2i+i/A2i+2(x^*'''^) = det [Aii[X^ , ■ ■ ■ , x~ 



2n-l^ 

Ai ... A2"-i j=0 



= nn(i-<^'^)"'- (39) 

i>j a>b 

5.3 Skew Schur function representation of generalized scalar prod- 
ucts 

Recall that matrix elements of the monodromy operators are expressed in terms of the 
skew Schur functions (fT7|) and (fTHj). From Proposition 15 -H the generalized scalar product 
of the form (1231) is written as 



n (ttU)^ n \uJ)M E ((^^")-^)^AVA<(^^")^) ■ ■ ■ 

j=l Vll ^ fc=i Vll ^ A'i<A/, l<i<2{n+m)-l 

X SA2n-l/A2n-2((M'"^)~^)sA2n-l/A2n((M'^^)^) X S A2n + 1 / A2" ( ("U^^ ) ~^ ) S A2n + 1 / A2n + 2 ( (m"^^ ) ^ ) 
X Sx2{n+m.)-l /A2(™+m)-2((n"^'")~^)SA2{n + m)-l((n'^'")^) 



where (u'^")"'^ = ((u^")"^, ■ ■ ■ , and so on. In particular, consider the plane 

partition in an x L x oo -box. This corresponds to (0| ULl-n+i C{q-^'^) Ilfcf A{q-^l^) 
HfcLi -S(q'*-'^"^''^^)|0). If we substitute the parameterization (jHH) into the above relation, we 
may reproduce with Af — oo by the help of (jHHl), as expected. 



6 Direct Calculation of Generalized Scalar Products 

In the previous section, we have established a determinant formula for the generalized 
scalar product in terms of the skew Schur functions. In this section, we give a more 
explicit determinant formula. This is done by a direct application of the relations obtained 
from the QISM in Section 12.11 The advantage of the formula here is a simplicity for 
evaluations. Actually in the next section, we demonstrate an evaluation of the determinant 
corresponding to the boxed plane partition in an x L x M-box. 

We remark about the notations. In this section, we often drop the indices of vari- 
ables to avoid notational complexity if there is no confusion. For instance, we abbreviate 
Y^=iK{'^t) ^'^ ^{'^^) '^^ j^s^ ri-^'l^"^) "^^^ Define the functions A(u,w) = 

m2 - v\ a{u) = and d{u) = u^'+\ 
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Theorem 6.1. Generalized scalar products have the following expression in terms of a 
determinant. 



N„ 



Ml 



(oi n cinf") n ^(^^) ■ ■ ■ n ^'(^f ^) n 

i=l i=l i=l 1=1 

Li Kl Lm 

n ^) \{B{uf^) w D{uf-) n B{uf-m 

i=l i=l i=l i=l 

i<j i>j 

(-l)^^Met G ({m^^^}; {m^^""}) (40) 
w/iere G = (c^j^ is an (N + M + L) x (N + M + L) matrix with 



a = 0, 6 = 0, 

l<i<A^ + M, 1<J<K + L{= N + L) 

a = 0, 1 < b < n, 

1 <t < N + M, 1 < J < Mb 

1 < a < m, 6 = 0, 

l<i< La, 1 <3 < K + L 

1 < a < m, 1 < b < n, 



(41) 



and the functions are 



h,^^{u) = J] [J] A(n^^n)A(n^^n)] a{u)u~'^^-'^~'^'^'^ 



(42) 



/^a,.(«)=nin^("'"''')^( 

i=l 



u, u 



.D, 



d{u)u 



2{i-l)+E"laKe-EZt^- 



and 



n (n 

i=l L 



I I I ^A,, \ 



n 



E}=iM,-E".,+iM,-l 



(43) 



n (n 

1=1 L 



n 
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Here we set N = Y^^i, M = Y^^i, K = Y^Ki = N and L = ^L,. The se- 
quence {ut''^}lY' is arranged as {uf"}.="i followed by {ut"}fi\, {uf 
Similarly, the sequence {wf^*"}^!^ is arranged as {wf'"}^™ followed by {wf""}^]^, ■ ■ ■ , 

Proof. Firstly, since © holds, the generalized scalar product is symmetric in each of {u^'}, 
{u^'}, {m*"'} and {u^'}, respectively. We can generalize the relations (|Tn|) and (fT!?jl to 

AT M 

= E(n^(-^''-^'))" (n-^') " (n-^)n^(-^')n^(-^')- (44) 

A',C' 

Here the summation is taken over all decompositions {uf'}^i U {uf'}fLi = {ufjfL^ U 
{uf}f^^. The relation pijl is obtained as follows. First, the coefficient of Y[ A{u^) 11 
is clearly Y[ f{u^i u^)- To see the coefficient of Y[ A{u'^') Y[C{u^'), exchange the variables 
of nC'(^i'^) n^(^^) by relation ^ till H 11 ^(^^')' then move the A's through 

C's and we have a coefficient 

V fiuf'^.K^)^^ 

where each uf^' E {u'"} fl {u^'} and uf-^' G {W^} H {u'^'} appears once in the product. By 
the definition (jSI) of / and g, we arrive at (jUj). 

Fix the number of the operators C and B both A^. For the case where there is no 
y4's and no -D's, the formula is nothing but the scalar product (|2T| . Assume the formula 
for n and m, that is, n products of A's and m products of D's are inserted. Then we 
calculate the case of n + 1 insertion of products of A's. Here we set A^' = + ■ ■ • + A^„, 
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M' = Ml + ■ ■ ■ + M„. From ()44|) and the formula for n and m, we get 



(0| Y[ C{u^"+') Y[ A{u^"+')---Y[C{u^')Y[A{i 

Ll Kl Lm Km 



,N 2A/„+i-l 



A',C' 

M„+i N„+i N„ M„ Ni Ah 

(oi n ^(^^') n c(n^')nc'(^''")n^(^''")---nc'(^''^)n^(^^^) 

Ll Kl Lm Km 

Yl Diu""') Yl ■ ■ ■ n Diu"""^) Yl 5(n^'")|0) 

i<j i>j 

detG({M^"}; JJ «("^') ^ ^ (45) 

where {wf'j^'J^^^ = {uf'^'~^}fj[^^\{uf}^^^\ and the summation is taken over all decom- 
positions {m^'} U {u'-^'} = U Here we used {0\A{u) = {0\a{u). We now 
calculate X, 

1 /„ „a2M„+i-1 /„ ^ 



i>j \i>j 



i<j \ i<j 

X 4 (Kn, ■ ■ ■ , {u^"} u K'}; Kn, ■ ■ ■ , K"}) 

x'.n^ ) (n^ 



Cn + l 
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We deal with the last expression in detail. First, 

i<j \ i<j / 



-ly n A(«f , o n in «^') n ^(^^-^ 

i<j j=l 

N Ar'+7V„+i+2(Af'+M„+i)-2 



■ir 11"' 



det 



i=l 



-(7V'+iV„+i+2M'+2(^-l)) 



l<k,l<Mn + l 



In the last line, we used the Vandermonde determinant. Here, a reorders {u^'^'~^} into the 
sequence {u'^'] followed by {u"^"}, {m"^"}, ■ ■ ■ , {u^^}, {u^'}- Next, 



n«-0 (n 



7V'+Af„+i+A/'-l 



,N Ar'+Af„+i+2(M'+M„+i)-2 



(Ar'+Af„+i+Af'+Af„+i-l) . 



U 



c 



2A/„+,-l 



E}.iM,-E^.+iM,-i 



c 



M' 



A^'+JV„+i+A/'~l 



u 



c 



n 

{Ar'+Ar„+i+A.f'+Af„+i-l) 



:n 



u 



2A/„+i-l 



,N Ar'+Ar„+i+2{A/'+Af„+i)-2 
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X 



n (H " 

i=l L 

n 



u 



-M: 



n + 1 



,N A/'+A/„+i-l 



n 



Cn + l 



i=l 



n (H "" 



n 



X n 



,C„+i 



(n 



In + l 



M' + M„+l-l 



n+1 r 



n 



Note that H'"'^'!!'"*^' = H ""^""^^ 11 '"'"^""^^ • Going back to substituting the above 

calculations, we get 



i<j 



({«^}; {u^}) i-ir'' (-l)-det G [{u^y, {u^-y) H 

.,N ~{N'+N„+i+2M'+2{i~l)) 



A',C' 



det 



1=1 



l<k,e<Mn+l 



i<j i>j 



A',C" 



]V+]\/' + l<f<]V+A/'+A/„+l 



where a' sends the indices of row with {n^'} for matrix G ({u^^^}; {n-^^^}) to the last. 
There appears another sign (— because L lines of the matrix G with 1 < a < m 
are passed through by M„+i lines. We can extend the summation of A', C to all {w^'} C 

becomes 



= unless {u"^'} C U {u'^"+^}. As a result, the sum 
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The whole output is exactly the n + 1 case of the formula. Computation for m + 1 case is 
done in the same way by developing relations and (|THl . Then by induction the proof 
is completed. □ 



7 Generating Function of Boxed Plane Partition 

As an application of the formula in the previous section, we evaluate the determinant of 
the generalized scalar product corresponding to the generating function of the boxed plane 
partition in an x L x M-box. 

Proposition 7.1. The following product expression holds for the generalized scalar product. 

L L-N N 

(0| n C{q~'l')\{A{q-''^)\{B{q^'~'y'm 



k=L-N+l k=l k=l 

^ ^-lLMN+^iL~N){M+l){L~N+l) ^ ■ ■ i (46) 



i=l j=l 



Proof. For an evaluation of the relevant determinant in ()4U|). we use the following deter- 
minant formula extending the one in fH] : 



det 



^(i+j-l)/2 _ ^-(i+j-l)/2 

g(i+i-l)/2 _ g~{^+j-l)/2 J ^^^^^ yy--^ > ) 1<^^L 



l<j<N 



l<j<L~N 



^_^-^lN{N-l)gliL-N)iL^+N^-LN+L-N)^^(L-N)(L-N+l)^gl/2 _ ^-1/2^) iL(L-l)+|7V(JV-l) 
N-1 L-N L-1 

Hi- - kf-' Hi- + n 



X 

k=0 k=l k=L-N+l 



X 



n [^-^1' n [^-^1 n [^+3-^" (47) 

l<1<i<Af Af+l<7<j<L L-N+\<i<L 

-J - l<j<JV 

where we put s = q'^ = and [k] = ^1/21^-1/2 • The remaining computation for the 

proof is a straightforward exercise. 

Eq. ()47|) is briefly proved as follows. The determinant of the matrix multiplied by s*/^ 
for rows with 1 < i < L and s^-?"^)/^ for columns with 1 < j < is a polynomial 
in s of degree LN + {L - N){L - N + l)/2. If we put s = g^, < A; < A^ - 1, the 
matrix is spanned by linearly independent vectors Xj = (^g(2i+i-fc)i/2j ^ I < j < k, and 

~ (^g(2-?+^+^)*/2) , I l£ j 1^ L — N, and its rank is k + L — N, so we have the divisor 
{s—q'')^~^. If we put s = q~^, I < k < L—1, the matrix is spanned by linearly independent 
vectors Xj = (g(2i+i-fc)«/2~j ^ 1 < j < max{/c, L — A^}, and its rank is max{k, L — N}, so we 
have the divisor (s — for 1 < A; < L — A^ and {s — q~^)^~^ ioiL — N+l<k<L — l. 
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We only need the overall factor. It is obtained by examining the coefficient of the lowest 
degree term of s. It turns out to be the product of a Cauchy-type determinant and a 
Vandermonde-type determinant, so it can be easily calculated. □ 

8 Conclusion 

In this paper, we have introduced and calculated a generalized scalar product of the inte- 
grable phase model. We have shown in two ways that a generalized scalar product with 
generic spectral parameters is expressed as a determinant. First, we have shown that the 
matrix elements of the monodromy operators are expressed in terms of the skew Schur 
functions through graphical interpretations; vertical lattices with arrows and lattice paths 
on them. A generalized scalar product is expressed as the sum of products of the skew 
Schur functions, where the sum is taken over all the partitions satisfying a restriction 
Ai < M. The sum is reexpressed as a determinant, applying the Cauchy-Binet formula. 
Second, the commutation relations for the monodromy operators of the QISM allow us to 
express a generalized scalar product as a determinant of matrix with elements of two-point 
functions and one-point functions. This expression is also valid for the inhomogeneous 
phase model. 

We have found that a generalized scalar product is related to the generating function 
of the skew plane partition. More precisely, a generalized scalar product with the special 
choice of the spectral parameters is up to some factor equal to the generating function of the 
corresponding skew plane partition. A sequence of the monodromy operators determines 
the skew shape. 

The skew plane partition in the semi-infinite box was considered in terms of the free 
fermion representation in [3 E] • The generating function of the plane partition can be writ- 
ten as the expectation value of products of successive vertex operators. The phase model 
representation here of the skew plane partition has some advantages when restricted in a 
box. In the fermion formalism this restriction requires us to insert projection operators 
between vertex operators and makes it difficult to calculate the generating function explic- 
itly. Correlation functions are easily treated in the fermion formalism. To find an explicit 
expression of general correlation functions in the phase model formalism is remained as a 
future problem. 
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